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Abstract. 

The stability of an attractive Bose-Einstein condensate on a joint one-dimensional 
optical lattice and an axially-symmetric harmonic trap is studied using the numerical 
solution of the time-dependent mean-field Gross-Pitaevskii equation and the critical 
number of atoms for a stable condensate is calculated. We also calculate this critical 
number of atoms in a double-well potential which is always greater than that in an 
axially-symmetric harmonic trap. The critical number of atoms in an optical trap 
can be made smaller or larger than the corresponding number in the absence of the 
optical trap by moving a node of the optical lattice potential along the axial direction 
of the harmonic trap. This variation of the critical number of atoms can be observed 
experimentally and compared with the present calculation. 
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The successful detection fT]!^ of a Bose-Einstein condensate (BEC) of dilute weakly- 
interacting trapped bosonic atoms at ultra-low temperatures initiated intense theoretical 
activities on different aspects of the condensate. The experimental magnetic trap is 
usually axially symmetric. More recently in different experiments a periodic optical- 
lattice potential generated by a standing-wave laser field has been employed along the 
axial direction of the magnetic trap jSHHEUHl HI IE] The optical-lattice potential has 
been used in the study of matter- wave interference [3], of oscillating atomic current in 
a one dimensional array of Josephson junctions [Sj, of Bloch oscillation and Landau 
Zenner tunneling and of superfluid-insulator classical [7] and quantum [I] phase 
transitions among others. The periodic optical-lattice potential has also been useful in 
the generation of matter- wave bright soliton jH] . There are many theoretical studies on 
a BEC in a periodic optical-lattice potential [Oj. 

Most of the above experiments were performed employing atoms with repulsive 
interaction, where the BEC is stable for any number of atoms. However, atoms with 
attractive interaction have been used in some experiments [21IBIEIIIE] , where the BEC 
is stable for the number of atoms smaller than a critical number When 
the number of atoms increases beyond this critical value, due to interatomic attraction 
the radius of the BEC tends to zero and the central density of the condensate tends to 
infinity. Consequently, the condensate collapses while emitting atoms in an exploding 
fashion due to three-body recombination until the number of atoms is reduced below the 
critical number and a stable configuration is reached [2] . With a supply of atoms from an 
external source the condensate can grow again and a series of collapse and explosion can 
take place and has been observed in a BEC of 7 Li atoms [2]. The attractive condensate 
has been fundamental in the generation of matter-wave bright soliton [8 j and in the 
study of collapse and explosion of a condensate simulating the supernova explosion of a 
star [To] . 

The critical number of atoms in an attractive condensate in an axially-symmetric 
trap has been studied experimentally [311011111 and numerically (T^HUE] by various 
authors under different symmetries extending from a pancake shaped BEC through 
spherical to a cigar shaped one. When the axial trapping frequency ui z equals the radial 
trapping frequency uj p one has a spherical BEC. When u z » u p (<< ui p ) one has a 
pancake (cigar) shaped BEC. In view of recent experiments [SI EJ El El HI IS] on a cigar- 
shaped BEC trapped jointly in a axially-symmetric harmonic plus a optical-lattice trap, 
we investigate the stability of an attractive cigar-shaped condensate in such a trap using 
the direct numerical solution of the mean-field Gross-Pitaevskii (GP) equation [16J with 
appropriate trapping potential without any simplifying assumption. 

Specifically, we calculate the critical number of atoms in an attractive BEC trapped 
jointly in an axial and optical-lattice traps. This will have direct consequence on the 
study of the generation and movement of bright solitons as in the recent experiment [B] 
with similar traps. We shall see that this critical number could be increased by 
appropriately setting the optical-lattice trap in the axial direction. A larger critical 
number is always welcome from the experimental consideration. Also, as the results 
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of the present mean-field study could be verified experimentally, this will provide a 
stringent test of the applicability of the GP equation on the study of collapse of an 
attractive condensate in the presence of a joint optical plus axially-symmetric traps. 
The most recent experiment JT] on the critical number in an axially-symmetric trap 
alone agrees well with the theoretical predictions [T3*|IT1]. 

There have been many theoretical studies of a BEC in a double-well potential 
and an experimental study of the same seems to be under control ^H]- in view of this 
we also present a study of the critical number of atoms in an attractive condensate in a 
double-well potential trap. 

The time-dependent BEC wave function ^(r; r) at position r and time r is described 
by the following mean-field nonlinear GP equation ^H] 



d h 2 V 2 

"^~^T + vr(r) + ^ iV|vI/(r;r) 



*(r;r) = 0, (1) 



where m is the mass and N the number of atoms in the condensate, g = 47rh 2 a/m the 
strength of interatomic interaction, with a the atomic scattering length. In the presence 
of the combined axially-symmetric and optical-lattice traps V(r) = \muj 2 (p 2 + v 2 z 2 ) + 
Vaux where u is the angular frequency of the harmonic trap in the radial direction p, vuj 
that in the axial direction z, with v the aspect ratio, and Vaux is an auxiliar potential 
simulating a double-well or optical-lattice trap introduced later. The normalization 
condition is J dr|\E'(r; r)| 2 = 1. 

In the axially-symmetric configuration, the wave function can be written as 
\&(r, t) = if)(p,z,r). Now transforming to dimensionless variables x = y/2p/l, y = 
y/2z/l, t = tuj, with the harmonic oscillator length I = ^Jhj (mui), and ip(x,y;t) = 

xJ¥J\/i^{p,z]T), ((H) becomes [13] 

d d 2 ld_ d 2 1 ( 2 2 2 v 



dt dx 2 x dx dy 2 4 



Vaux 1 „ rx 

+ — ^ + 8V2nn 

nuj x z 



<p(x,y;t) 



2 



<p(x,y;t) = 0, (2) 

where nonlinearity n = Na/l. In terms of the one-dimensional probability P(y,t) = 2n 
J °° dx\<f(x, y, t)\ 2 /x, the normalization of the wave function is given by dyP(y, t) = 
1. 

A double-well potential for an axially-symmetric trap can be simulated by taking 
Vaux _ V^ w 2 

where V is the strength of the double- well potential V^ w . The optical potential 
created with a standing-wave laser field of wavelength A along the axial direction can 
be represented by taking Vaux = V^p^ = V E R cos 2 (k L z), with E R = h 2 k 2 L / '(2m), 
k L = 2tt/X, and V the strength jHllllElin]- This periodic potential has a maximum 
at the center of the harmonic trap z = 0. The center of the harmonic trap could 
also be made to coincide with a minimum of the periodic potential or anywhere in 
between. All these cases can be covered by taking a phase 5 in the cos 2 (kiz) term 
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Figure 1. The reduced critical number n cr ^ vs. strength Vq of the double- well 
potential for v = 0.5: line with H > double well; straight line n CT ^ — 0.63 — > Vq = 0. 



so that V r 



opt 



VqEr cos 2 (kiz + 5). In terms of the dimensionless laser wave length 



A = y/2X/l and a dimensionless standing-wave energy parameter E R /(ftw) = 47t 2 /Aq, 



We shall see that by appropriately choosing the value of 5 and Ao the critical number 
of atoms in an attractive BEC could be increased after the introduction of the optical 
trap. 

We solve (J2J numerically using a split-step time-iteration method with the Crank- 
Nicholson discretization scheme described recently [T^j. The time iteration is started 
with the harmonic oscillator solution of (J2J) with n = 0: p(x,y) = [z^/ (87r 3 )] 
xe -{x 2 +vy 2 )/i ^| TYie nonlinearity n and the optical-lattice potential parameter Vq 
are slowly changed by equal amounts in 5000n steps of time iteration until the desired 
value of n and Vq are attained for the double-well or the optical-lattice potential. Then, 
without changing any parameter, the solution so obtained is iterated 5000 times so that 
a stable solution is obtained independent of the initial input and time and space steps. 

For a repulsive condensate the scattering length a is positive. For an attractive 
condensate the scattering length a and the nonlinearity n are negative. In the calculation 
we find that, for an attractive condensate with negative n, stable solution of the 
GP equation cannot be obtained for |n| greater than a critical value ^ cr j^ — the 
reduced critical number. For the spherically symmetric 0.5746 
The critical value for the number of atoms is given by A^ cr ^ = n CT ^l/\a\. The 
dimensionless parameter n QT ^ is universal, whereas the number iV cr j^ depends on 
the typical experimental set up, e.g., the atomic scattering length and the harmonic 
oscillator length of the trap. 



the periodic optical-lattice potential V^pt of (J2I) is 




(4) 



Critical number of atoms in an attractive condensate 



5 




Figure 2. The reduced critical number n CI u vs. dimensionless wave length Ao for 
v = 0.5 and (a) V = 1, (b) V = 4, (c) V = 8, and (d) V = 12: line with + ^5 = 0; 
line with * — > S — tt/2; straight line n cr ^ — 0.63 —> Vq = 0. 

First we present a study of the critical number n CT u for the double-well potential. 
We consider an elongated cigar-shaped harmonic trap with v = 0.5 although the 
calculation can be extended for other values of v. In figure 1 we plot n CT ^ vs. Vq. 
The case V = corresponds to the absence of the optical-potential or double-well trap 
and in this case we obtain n cr j^ = 0.63 in agreement with k cr ^ = n^^u 1 ^ 6 = 0.56 
reported in reference ^3] for v = 0.5. We find from figure 1 that ^ cr j^ for the double- 
well potential is always larger than that for the corresponding harmonic axial trap. The 
double-well trap separates the original condensate effectively into two parts with a region 
of low density in between. When Vq is large this separation is more pronounced and 
one effectively has two separated condensates and the critical number of atoms in the 
double-well trap tends approximately to twice the critical number in the corresponding 
harmonic trap alone. For V = the double-well is absent and ra cr j+ = 0.63. Figure 
1 shows the evolution of ^ cr jt between these two limiting values. It should be recalled 
that the exact limiting value for n cr ^ for large Vq as well as the evolution between 
the two limits will depend on the type of the double-well potential employed. Figure 1 
represents these general features for the double- well potential ©• 

In most experimental set ups with an optical-lattice potential a cigar shaped trap 
has been used and we consider only this symmetry [y < 1) in this paper. Here we 
consider the anisotropy parameter v = 0.5 and four values of the strength Vq of the 
optical-lattice potential: Vq = 1,4,8 and 12. We consider the two limiting cases 5 = 
and tc/2 in the optical-lattice periodic potential (jlj). 

The critical nonlinearity n CT ^ vs. dimensionless wave length Ao for different Vo 
and 5 (= 0, 7r/2) are plotted in figures 2. This constant critical value (n cr j^ = 0.63) 
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Figure 3. The full axial potential V ax j vs. y for Vo = 8, 6 — and v — 0.5 for (a) 
A = 2, (b) A = 8, (c) A = 12, and (d) A = 20: full line -> y axi ; full line with x -> 
axial harmonic potential v y 



for z/ = 0.5 and Vo = is also shown in figures 2 for a comparison. The lines with 
+ represent 5 = or the cos 2 (27n//A ) dependence of V p^ whereas the lines with * 
represent 5 = n/2 or the sin 2 (27ry/Ao) dependence. For a fixed Vo and 5 = one 
has n. cr j+ < 0.63 for small A . With the increase of A , ^ cr jt increases past 0.63 and 
eventually it decreases to 0.63 asymptotically for a large enough A . For 5 = n/2 and 
small Ao, ^ cr jt starts at a value very similar to the 5 = case. With the increase of 
A , ^ cr it monotonically increases in this case to the value 0.63, asymptotically, however, 
remaining always smaller than 0.63. The plots of n CT ^ vs. Ao for any 5 lie between the 
5 = and 5 = tc/2 limits shown in figures 2. 

The behavior of the n CT ^ vs. Ao curves reported in figures 2 can be understood 
from a consideration of the total axial potential 



^axi = \" 2 y 2 + Vo^r 



cos ( — y + 5 



(5) 



It should be noted that in the presence of the periodic optical-lattice potential the radial 
trap remains unchanged as the optical lattice does not provide any confinement in the 
radial direction. In figures 3 we plot V ax j vs. y for four different values of wavelength 
Ao for Vo = 8 and 5 = 0. For Ao = 2 we have a large number of nodes in the potential 
V ax j revealing a large number of confining traps of the periodic potential as one can 
find in figure 3 (a). The number of available traps of V ax j in a certain region of space 
in the axial y direction decreases as Ao increases (to 12 through 8) as found in figures 3 
(b) and (c). Eventually, for a large enough Ao (~ 20) the confining traps of the periodic 
potential V ax j disappear and it reduces essentially to the harmonic trap potential as in 
figure 3 (d). For large Aq the extention of each of the periodic traps of V ax j is larger 
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Figure 4. Same as in figures 3 but for 5 = tt/2 

than that of the axial harmonic trap, and V^ x j tends towards the axial harmonic trap 
u 2 y 2 /4. Hence one can safely conclude that for large A , ^ cr jt tends to its value (= 0.63) 
in the absence of the optical potential as one can find from figures 2. 

The typical axial dimension of the attractive condensate before collapse is larger 
than the width of a single confining trap of the optical potential for small Ao (figure 3 
(a)), hence the condensate should occupy a few of these confining optical-lattice traps. 
Consequently, the BEC will be squezeed towards the center of each these confining traps 
with a region of low density between the traps. Because of this further squeezing the 
BEC will collapse with a smaller number of atoms and lead to a n CT ^ smaller than 
that in the absence of the optical trap. This effect will be large for a stronger optical 
potential with large Vq as one can see in figures 2. The typical dimension of the attractive 
condensate in the axial direction is smaller or comparable to the dimension of a single 
trap of the optical potential for larger Ao- For 5 = the appearance of a maximum of 
the axial potential V^ x j at y = pushes the condensate apart into two nearby minima of 
the axial potential. Consequently, in these cases the BEC swells and one effectively has 
two separated pieces of the condensate in two traps. In this case the central part of the 
total trap, where the BEC is formed, is similar to a double-well potential. Consequently, 
there is more space available for the full condensate and hence n CT ^ can grow beyond 
the critical value in the absence of the optical trap. 

In figures 4 we plot V ax ^ vs. y for four different values of wavelength Ao for Vq = 8 
and 5 = tt/2. In this case one has a minimum of the axial potential V^ x j at y = in 
contrast to the 5 = case where one has a maximum of at y = 0. Consequently, for 
larger Ao the BEC will be essentially squeezed into the central trap of the optical-lattice 
potential. The central part of the total potential, where the BEC is formed, is never of 
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the double-well type, but rather like a stronger single-well type. Because of the above 
squeezing the critical number of atoms will always be less than the critical number in 
the absence of the optical-lattice trap. Again n cr ^ grows monotonically between the 
small and large Ao limits as one can see in figures 2. 

In conclusion, using the numerical solution of the GP equation we have studied the 
stability of an attractive BEC trapped jointly by an optical and an axially-symmetric 
harmonic traps. We calculate the critical number of atoms for stability in this case. 
We find that the reduced critical number n cr j^ = N CT ^\a\/l could either increase or 
decrease in relation to the same number in the absence of the optical trap. However, 
for most experimental set ups Ao is expected to be less than 3 and and from figures 2 
we find that n CT ^ should reduce after the introduction of the optical lattice trap. For a 
double-well potential n CT ^ is always larger than the same for the corresponding axially- 
symmetric harmonic trap. A prior knowledge of the critical number n cr j^ as obtained in 
this paper will be useful in the experimental generation and study of matter-wave bright 
solitons under similar trapping conditions. Also, n CY ^ can be measured experimentally 
and compared with the present theoretical results. This will be helpful in evaluating 
the applicability of the mean-field GP equation in the study of attractive condensates 
trapped in a joint periodic optical plus a harmonic potential. 
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